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ABSTRACT 



A fraction of the radiation produced by an accretion disk may be Thomson scat- 



OO 

, tered by a wind flowing away from the disk. Employing a simple plane-parallel model of 

the wind, we calculate the polarization of the scattered radiation and find that its sign 
I depends on the wind velocity, (5 = v/c. In the case, 0.12 < /3 < 0.78, the polarization 

' is parallel to the disk normal, i.e., it is orthogonal to the standard Chandrasekhar's 

polarization expected from accretion disks. The velocity of an wind is likely to sat- 
urate near the equilibrium value /?^, ~ 0.5 for which the accelerating radiation pressure 
is balanced by the Compton drag. Then the change of polarization by the wind is most 
\ pronounced. This may help to reconcile the standard accretion disk model with the 

fSJ ! optical polarimetric observations of non-blazar AGNs. 
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' 1. Introduction 

The standard model of a black hole accretion disk predicts the polarization of the emerging 
radiation. The radiation gets polarized due to electron scattering that strongly dominates the 
^ ' absorption opacity in luminous accretion disks. Polarization of light emerging from an electron 

d ' plane-parallel atmosphere was calculated by Chandrasekhar (1960) and Sobolev (1963). According 

to their results, the polarization should be parallel to the disk surface with a maximum of 11.7 % 
when the disk is viewed edge-on. This simple model, however, disagrees with observations: the 
optical polarization in non-blazar AGNs is typically ~ 1 % or less, and it tends to be parallel to the 
radio jet that is presumably perpendicular to the disk (Stockman, Moore, & Angel 1984; Antonucci 
1992). The problem of accretion-disk polarization has been discussed by a number of authors (e.g. 
Gnedin & Silant'ev 1978; Loskutov &: Sobolev 1981; Sunyaev & Titarchuk 1985; Phillips & Meszaros 
1986; Coleman & Shields 1990; Laor, Netzer, & Piran 1990; Matt, Fabian, & Ross 1993; Kartje 
1995; Agol & Blaes 1996). In this Letter, we suggest that the original polarization of the disk can 
be changed due to Thomson scattering in a mildly relativistic wind. 
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Gas outflows axe observed in non-blazar AGNs, and evidence for mildly relativistic bulk veloc- 
ities has been reported (Lcighly et al. 1997). The outflows are seen on scales far exceeding the scale 
of the central source ^ 10^'^ cm, and the gas dynamics in the central region is poorly understood. It 
is quite probable that there is a wind from the innermost region of an accretion disk where the bulk 
of the observed radiation originates. The wind may form due to a gas outflow from a corona of the 
disk. It may also form due to an outflow of electron-positron pairs produced by the gamma-rays 
from the disk. In the latter case, the created pairs cool efficiently by inverse Compton scattering 
down to the Compton temperature kTc ~ 1 — 10 keV, and the cool light plasma is then easily 
pushed away by the radiation pressure. 

To investigate the effects that a wind could introduce in the pattern of polarization, we consider 
the simplest plane-parallel model in which the disk is replaced by a homogeneously emitting plane 
and the wind is assumed to be a cold vertical outflow from the plane (by "cold" we mean that 
the thermal motions of particles in the outflow are slow compared to the bulk motion). Thomson 
scattering of frequency-integrated polarized radiation is governed by the transfer equations written 
down in § 2. The equations account for the relativistic aberration of light in the wind rest frame, 
which is crucial for the polarization due to scattering. In § 3 and § 4 we discuss in turn the optically 
thin and optically thick winds. The results are summarized in § 5. 



2. Basic equations 

In a plane-parallel slab, the polarized radiation is represented by the frequency-integrated 
intensities Ii,Ir, where the index I refers to the radiation polarized in the meridional plane (defined 
by the normal and the ray) , and index r refers to the polarization perpendicular to this plane. The 
radiation is axisymmetric with respect to the normal, with the intensity being a function of /x - 
the cosine of the angle between the ray and the normal. Stationary Thomson transfer in a cold 
medium moving in the vertical direction with a velocity f3 = v/c is described by the equations, 

(97 do 

= (1 - /?//)(/ - s), = (1 - P^^){Q - R)- (1) 

Here, r is the scattering optical depth that the medium would have at rest, / = Ir + lu and 
Q = Ir — h- S and R are source functions representing the scattered radiation. 

A simple way to derive the source functions is to transform the radiation field to the wind 

rest frame where Thomson scattering is coherent. Polarization is invariant with respect to Lorentz 
boosts along the normal, and the frequency-integrated intensities I and Q transform as (see, e.g., 
Rybicki & Lightman 1979) 
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where the index c stands for the comoving frame, D = 7(1 — (3fx) is the Doppler factor, 7 is the 
Lorentz factor. When viewed from the comoving frame, the scattered radiation is represented by 
the source functions, S'^{p,c) and R'^^ic), given by Sobolev (1963). Transforming the scattered 
radiation back to the laboratory frame, S{fj,) = D~*S'^{iic) and R{ii) = D~'^R'^{jj,c), we get 



^(/^) = - lis + Q'o- Q^), 

where I^, ("^ ~ 0' ^) moments of the radiation field in the comoving frame. 



(2) 
(3) 
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are related to the corresponding moments in the lab frame, 7^, by 
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-/3(/o + /2) + (l + /3')/i 



Qot transform in exactly the same way. Note that — Q2 = Qo — Q2- 

In a 2D approach, an initially vertical wind becomes quasi-radial at some height and its opacity 
falls off. We consider the simplified ID problem in which the wind is replaced by a slab of vertically 
outflowing plasma of a given optical depth tq. At the outer boundary, r = 0, we assume free escape 
of the radiation. At the inner boundary, r = tq, we assume the existence of a source of limb- 
darkened positively polarized radiation (Chandrasekhar 1960) for which Ii « 0.581Io, I2 ~ 0.414/o- 



3. Optically thin wind 

In the case of a small optical depth, tq, the optically thin approximation holds at angles 
satisfying the condition, t^j, = tq{1 — (5^)/ iJi < 1, and equations (1) yield the scattered radiation 
emerging from the slab as 

Isc ^T^S {jl) , Qsc'==' T/^RilJ-)- 

The polarization degree of the scattered radiation equals Pscil^') = Qsc/hc = R/S. Prom equations 
(2) and (3) we get 

^ 3(1-Mg)r? 

(8/3) + (3/x2-l)^' W 

where we have introduced the parameter 

^ _ ^|-(^73) + Qo-Q2 ^ 1 
-'0 -'0 ^ 
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The contribution to rj from Qo and Q2 is ~ 30 times smaller than that from I2 — Iq/3, and can be 
neglected. The parameter —1/3 < ry < 2/3 shows to what extent the radiation field is stretched 
along the normal in the rest frame of the wind, which is of crucial importance for the polarization 
of the scattered radiation as Psdl^) lias the same sign as r] at all In Figure 1 we show the 
dependence of 77 on the wind velocity (5: r] changes sign at Pi ^ 0.12 and P2 ~ 0.78. It follows that 
the scattered radiation remains positively polarized if the wind is slow, (3 < The change in sign 
of the polarization happens if the wind has a velocity in the range, f5\ < (3 < ^2- 

In the case the outflow is comprised of a light plasma, the time-scale to achieve Compton 
equilibrium with the radiation field is small and the wind velocity adjusts so that the radiation 
pressure is balanced by the Compton drag (see, e.g., Phinncy 1982). It means that, in the comoving 
frame, the net radiation flux vanishes, i.e., AttI^ = 0, and the equation for the equilibrium velocity, 

reads 

P.{h + h) = {l+[il)Iu (6) 

which yields /3* ~ 0.52 for the limb-darkened radiation from the disk. One can check that ?7(/3) has 
a minimum (5 = (5^. It follows that a wind produces the strongest negative polarization if it is 
in Compton equilibrium with the radiation field. This fact is illustrated by the diagram embedded 
in Figure 1 that shows the angular distribution of the disk radiation in the rest frame of the wind. 

A special feature of the aberration at mildly rclativistic (5 is that the radiation field is effectively 
"compressed" in the vertical direction. It is this compression that results in parallel polarization 
of the scattered radiation, and the compression is naturally strongest at /3 = for which the net 
flux in the comoving frame vanishes. 

Whatever velocity of the wind, (3, the scattered radiation has an extremum of polarization at 
fi = p. It corresponds to fic = 0, i.e., to the case the ray is perpendicular to the normal in the 
comoving frame, and we denote the extremum as p^. In Figure 1, is plotted against /3. For a 
mildly relativistic wind, r/, e.g., p^ = — 19 % for /3 = In the ultrarelativistic limit, the 

aberrated radiation gets concentrated in "head-on" direction and — 100 %, as discussed by 
Begelman & Sikora (1987). 

An observer viewing a disk with an optically thin outflow will see the sum of unscattered and 
once scattered radiation. In Figure 2 we plot the observed polarization as a function of the disk 
inclination, /i, for the case tq = 0.1 and /3 = = 0.52. The optical depth is small and the scattering 
in the wind weakly affects the intensity of the observed radiation. One can see, however, that the 
wind changes the pattern of polarization. This happens because the polarization extremum of the 
scattered radiation is boosted into angles fi ~ 0.5, where the original disk polarization is small. 
Note that the calculated slab model fails at small 12 as increases oc ij,^^ to = 1 where the single 
scattering approximation is not applicable. In a real wind, the scattering region has typical height 
comparable to its radius, and a detailed 2D model is needed to get exact pattern of polarization. 
To illustrate importance of geometry, take a cylinder of height h and radius r = V^h instead of 
slab. Then, at ;u > 0.5, changes to r^^ = \/3ro(l — /?m)(1 ~ A*^)"^^^- Combined with the slab 
source functions, S and R, this yields a rough estimate for the resulting polarization shown by the 
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dotted line in Figure 2. 

4. Optically thick wind 

An optically thick outflow may emerge in the case a fraction of the disk luminosity is 

emitted above 511 keV, as hard photons are absorbed by softer X-rays to produce pairs (the non- 
linear transfer of gamma-rays above a luminous disk is discussed in Beloborodov 1998). Compton 
equilibrium is established in each layer of the outflow, and the bulk velocity is determined by the 
radiation fleld itself through the equilibrium condition (6) . This means that the velocity should be 
calculated self-consistently in the transfer problem. Comparing the Compton equilibrium transfer 
to the classical transfer problem in an electron medium at rest, we note that in addition to a 
constant net flux, 47r/i(r) = const, we have l2iT) = const (one can check this combining the first 
moment of the transfer equation [1] with the equilibrium condition [6]). 47r/2/c equals the radiation 
pressure in the vertical direction, and its constancy reflects the condition that the net radiative 
force acting on vanishes in the limit /? — The equilibrium velocity varies with height, and 
pairs keep /? /3* as /?* is a strong attractor in phase space. 

We have simulated this transfer problem numerically. We take as initial conditions that the 
disk radiation is propagating freely in the slab, then "switch on" the scattering by the flowing pairs 
at t = 0, and follow the evolution of the system until a stationary solution is established. In the 

calculations, we use a grid which is homogeneous in the cos~^ fi and r directions. The number of grid 
points is Nj- x = 200 x 300. The stationary solution, /(r, n) and Q{t, fi), that we are looking for 
depends only on the parameter tq and does not depend on the density profile, n(r). Therefore, 
we choose the simplest profile, n^r) = const. The chosen time step equals At = O.bTo/naTcNr. 

The stationary solution for tq = 3 is shown in Figures 3 and 4. Figure 3 displays the wind 
velocity. In the lowest layers of the slab, the equilibrium velocity is less than 0.5c since besides the 
disk radiation there is some flux in the opposite direction due to radiation being back scattered by 
the outflow. In higher layers, the radiation gets collimated by the scattering on the moving pairs, 
and the equilibrium velocity increases to ^ 0.7c. Figure 4 shows the radiation emerging from the 
slab: it is signiflcantly anisotropic, and it is negatively polarized with an extremum oi p —14 % 
at « 0.43. 

5. Conclusions 

We conclude that due to Thomson scattering in a wind, the disk radiation originally polarized 
perpendicular to the normal can acquire parallel polarization. This effect arises due to relativistic 
aberration of light in the wind rest frame and takes place when the wind velocity is in the range 
0.12 < (3 < 0.78. Being characteristic for mildly relativistic winds, the production of parallel polar- 
ization does not occur in the previously discussed scattering by highly relativistic jets (Begelman 
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&; Sikora 1987) which produce strongly beamed radiation with polarization perpendicular to the 
jet. 

An observer viewing a disk through an optically thin wind will see the sum of the radiation from 
the disk and the radiation scattered in the wind. If the original disk radiation has the perpendicular 
Chandrasckhar-Sobolcv polarization then the scattering in the wind tends to diminish it. The 
strongest effect is produced by a wind with velocity /3 = /3* ~ 0.5, where /3* corresponds to a 
Compton equilibrium of the outflow with the disk radiation field. Such an equilibrium is likely to 
be established in a wind comprised of an plasma. We find in this case that a modest optical depth 
of the wind, ~ 0.1, is enough to change the pattern of the disk polarization: The polarization degree 
is reduced and there appears a range of inclination angles with negative (parallel) polarization. This 
change may naturally explain the optical polarization observed in non-blazar AGNs. 

We have also found that parallel polarization is produced in an optically thick outflow which 
might form above a black hole accretion disk. The outflow is semi-relativistic, and the emerging 
radiation is markedly beamed along the disk axis. The bulk of the disk radiation gets scattered, 
its original polarization is lost, and the emerging radiation acquires parallel polarization with an 
extremum of about —14 % at an inclination angle of ~ 7r/3. Note that an optically thick 6=*= wind 
also produces an annihilation feature above 511 keV, which has a blueshift up to ~ 2, albeit smeared 
out by the velocity gradient in the wind. 

Our general conclusion is that a wind produces parallel polarization when it is near Compton 
equilibrium with the radiation field. This conclusion is derived for vertical winds in slab geometry 
and it remains to be investigated how it is affected by deviations from this simple model. In 
particular, a magnetic field may aS"ect the wind direction, especially if magnetic pressure dominates 
the radiation pressure, and additional relativistic effects are important near Kerr black holes. A 
detailed 2D model should include a radial dependence of the disk emission and the wind density. 

I thank C.-I. Bjornsson, J. Poutanen, and R. Svensson for many useful remarks, and the referee 
for helpful comments. I acknowledge support from the Swedish Natural Science Research Council 
and RFFI grant 97-02-16975. 
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Fig. 1. — Dependence of parameter r] on the wind velocity, /?, (solid curve). The disk is modeled 

as a plane-parallel electron atmosphere emitting limb-darkcncd radiation (Chandrasekhar 1960). 
The dotted curve displays p^, the degree of polarization of the radiation scattered at right angle 
to the normal in the comoving frame. The embedded diagram shows the angular distribution of 
the aberrated radiation seen by an electron in the wind, I'^in)- Dashed, solid, and dotted lines 
correspond to /? = 0, /3 = /?* = 0.52, and (3 = ^2 = 0.78 respectively. 
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Fig. 2. — Angular dependence of the polarization from an accretion disk with a wind of optical 
depth tq = 0.1 and velocity = 0.52. The scattered radiation is calculated in the single scattering 
approximation which fails for the slab model (solid line) when the " observed" curve approaches the 
"scattered" one at ~ 0.1. The estimate for the cylinder model is shown by the dotted line. 
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Fig. 3. — Velocity profile of an e wind with optical depth tq = 3. 
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Fig. 4. — Radiation emerging from an optically thick e wind with tq = 3 in comparison to 
the Chandrasekhar's radiation of the disk: (a) - the angular dependence of the total intensity 
(normalized so that the net flux equals tt) , (b) - the angular dependence of the degree of polarization. 



